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ABSTRACT 


The  amount  of  Interaction  between  a  fluid  and  a  compliant 
coating  is  studied  for  a  one  dimensional  shear  fluid  disturb¬ 
ance.  A  thin  inhomogeneous  viscoelastic  layer  is  located  at 
the  interface  between  the  fluid  and  the  coating.  The  fluid  is 
assumed  to  have  no  mean  flow  field.  The  effect  of  different 
coating  properties,  thickness  of  compliant  coating  and  of 
transition  layer  as  well  as  the  frequency  of  oscillations  are 
analyzed . 
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I.  INTRODUCTION 


The  drag  of  a  body  can  be  decreased  by  a  compliant  coating 
only  if  a  considerable  amount  of  interaction  exists  between  the 
fluid  and  the  coating*  The  coating  has  to  be  able  to  be 
excited  by  a  fluid  disturbance  and,  in  turn,  the  coating 
disturbance  has  to  change  the  characteristics  of  the  fluid 
disturbance  in  such  manner  as  to  decrease  the  overall  drag  of 
the  body. 

In  typical  applications,  the  flow  field  will  have  either  a 
laminar  or  a  turbulent  boundary  layer.  For  low  Reynolds 
numbers  (air  application),  the  boundary  layer  will  typically  be 
laminar  whereas  for  high  Reynolds  numbers  (underwater  applica¬ 
tion)  the  boundary  layer  will  be  turbulent.  The  Interaction 
phenomenon  required  for  drag  reduction  may  be  totally  different 
in  the  two  cases.  The  drag  of  a  body  with  a  laminar  boundary 
layer  will  consist  mostly  of  the  drag  due  to  skin  friction. 

The  drag  of  the  same  body  but  with  a  turbulent  boundary  layer 
will,  however,  be  Influenced  the  most  by  the  magnitude  of  the 
Reynolds  turbulent  stresses.  The  drag  due  to  the  laminar 
sublayer  is  always  smaller  than  the  drag  due  to  the  normal 
Reynolds  stresses  acting  on  a  wavy  wall  (fig.  1). 

In  order  to  properly  model  the  fluid-structure  inter¬ 
action,  as  represented  by  the  flow  fields  discussed  above,  two 
types  of  fluid  disturbance  should  be  analyzed:  shear  disturb¬ 
ance  and  normal  (or  longitudinal)  disturbance.  The  shear 
disturbance  will  serve  to  simulate  the  Interaction  of  the  fluid 
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with  the  coating  for  either  the  laminar  boundary  layer  case  or 
the  laminar  sublayer  underneath  the  turbulent  boundary  layer. 
The  longitudinal  (or  acoustic)  disturbance  result  will  be  used 
to  model  the  interaction  between  the  Reynolds  turbulent 
stresses  and  the  compliant  surface. 

The  parameters  which  govern  the  local  properties  of  a  one 
dimensional  medium  can  be  taken  to  be  the  density  and  the  local 
wave  speed  of  the  medium.  A  disturbance  traveling  in  a  given 
direction  will  be  partially  reflected  whenever  a  discontinuity 
in  the  properties  of  the  medium  exist.  The  Interface  between  a 
coating  and  a  fluid  is  an  example  of  a  medium  discontinuity. 

If  the  effect  of  the  discontinuity  could  be  ameliorated  then 
the  reflection  and  the  transmission  of  energy  from  one  medium 
to  another  could  be  changed.  A  transition  layer  which  spreads 
out  the  discontinuity  between  the  two  media,  in  a  gradual 
manner,  can  be  expected  to  Increase  the  magnitude  of  the 
Interaction  between  the  two  media.  The  properties  of  the 
transition  layer  may  vary  continuously  from  that  of  a  solid  on 
one  side  to  that  of  a  fluid  on  the  other  side.  Alternatively, 
the  transition  layer  may  be  assumed  to  be  an  equivalent  layer 
whose  properties  are  determined  by  the  overall  local  motion  of 
the  surface  of  the  coating  (Fig. (2)).  The  transition  layer 
will  thus  exhibit  the  properties  of  an  Inhomogeneous  layer  with 
continuously  varying  properties  (Fig.  (3)).  Thus  the  transition 
layer  can  either  be  designed  as  part  of  the  compliant  coating 
or  it  can  be  thought  to  be  an  Integral  part  of  the  coating  with 
the  mean  transition  layer  properties  due  mainly  to  the 
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inability  to  proparly  identify  the  location  of  the  interface 
being  forced  in  a  random  manner  by  the  turbulence;  in  this 
context,  then  the  transition  layer  will  acquire  average 
properties  of  the  solid  and  the  fluid  and  the  actual  magnitude 
of  the  properties  will  depend  upon  the  average  of  the  space- 
time  history  of  the  surface  motion  at  each  location. 

In  the  two  previous  reports, 1,2  the  distribution  of  both 
shear  and  normal  disturbances,  due  to  a  disturbance  at  the 
bottom  of  the  compliant  coating  were  studied  (Fig.  (4)).  The 
results  indicated  that  the  transmission  of  shear  disturbances 
into  the  fluid  can  be  altered  considerably  by  the  presence  of 
thin  transition  layer.1  For  th®  same  conditions,  longitudinal 
(or  normal)  wall  disturbances  will  be  affected  by  the 
transition  layer  in  a  much  reduced  manner. 

Pedersen  et  al3  have  shown  that  an  inhomogeneous  layer 
with  an  exponentially  varying  Impedance  at  the  interface 
varying  from  that  of  a  solid  to  that  of  a  fluid  increases  the 
acoustic  energy  transmitted  from  the  solid  into  the  fluid. 
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II.  ANALYSIS 


The  enelysi8  consists  of  the  derivation  of  the  governing 
equations  for  the  fluid,  the  elastic  solid  and  the  transition 
layer  and  the  numerical  technique  utilized  to  integrate  the 
resulting  differential  equations. 

II. 1  Fluid 


Consider  a  one  dimensional  fluid  disturbance  located  at 
distance  h  from  the  lower  wall  (Fig  (3)).  The  equations 
governing  the  fluid  motion,  in  its  most  general  form  are  the 
Navier~S^okes  equations  and  the  continuity  equation. 
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where  ^  is  the  material  derivative  and  is  given  by 
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v.  and  v 0  are  the  local  fluid  velocities  in  the  x  and  the  y 
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direction  respectively  and  pf  i8  the  fluid  density.  The  fluid 
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where  y  is  the  fluid  viscosity. 

For  a  one  dimensional  fluid  shear  distrubance,  the  Vj 

component  of  velocity  term  will  be  zero.  If  the  fluid 
disturbance  is  assumed  to  vary  sinusoidally  In  time, 


v-j  *  v-j  exp(-iwt)  i 


then  the  governing  equation  reduces  to 


ay  (ua~ )  +  ipfwv1  =  o 


II. 2  Solid 


The  fluid  disturbance  will  propagate  from  its  source  to 
the  fluid-solid  interface,  where  the  no  slip  condition  and  the 
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continuity  of  stress  will  cause  the  disturbance  to  propagate 
into  the  fluid. 

The  generalized  equations  for  the  motion  of  the  solid  are 
the  Navler  equations: 

32U]  ^tHs  ^t12s 
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where  p$is  the  density  of  the  material.  The  solid  particle 
displacements  in  the  x  and  the  y  direction  are  u^  and  u2 

respectively.  The  solid  stresses  are  linearly  related  to  the 
particle  displacements  u^  and  U£» 
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where  G  is  the  shear  modulus  and  V  is  Poisson  ratio.  In 
general  G  will  be  a  real  number  if  the  solid  is  purely  elastic 
and  a  complex  number  if  the  material  has  damping  properties. 

If  the  fluid  disturbance  is  a  shear  vibration  and  is 
Independent  of  its  x  location  then  the  solid  will,  by 
necessity,  vibrate  in  a  similar  fashion;  thus  equations  10-14 


are  combined  into  an  ordinary  differential  equation 
h  dui 
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For  an  inhomogeneous  viscoelastic  medium,  the  shear 


modulus  G  will  be  a  complex  function  of  y. 


II. 3  Unified  Theory 


The  differential  equation  governing  the  fluid  motion  (Eq.  9) 
.and  the  equation  .  governing  the  solid  motion  (Eq.  15)  are 
seen  to  be  very  similar.  For  a  fluid  disturbance  with  no  mean 

flow  field,  the  fluid  disturbance  velocity  vi  *8  rel*ted  to  the 
particle  displacement  u^  by 

Vj  *  - iwUj  .  (16) 

The  vibration  of  the  fluid-solid  system  can  be  combined 
into  a  unified  equation  given  by 
A  dui 
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The  parameter  K  is  given  by 

G  solid 


K  - 


-iwy  fluid  , 


(18) 


and  the  density  p  is  either  the  solid  or  the  fluid  density. 

The  transition  layer  Interposed  between  the  fluid  and  the  solid 
will  be  assumed  to  vary  exponentially  from  the  shear  modulus  G 
on  one  side  to  the  fluid  viscosity  y  on  the  other.  The 


discontinuous  parameter  K  can  thus  be  written  as 


K ( y )  =  -ipw  +  [G  +  iyai]  exp(-  y/Jt)n  .  (19) 

For  values  of  y  close  to  the  inner  wall,  the  elastic  solid 
properties  are  recovered  and  for  values  of  y  far  removed  frees 
the  interface  region  (i.e.  y  >>£),  the  fluid  properties  are 
obtained.  In  equation  19,  i  is  the  nominal  thickness  of  the 
compliant  coating  and  n  is  Inversely  related  to  the  thickness 
of  the  transition  layer.  In  the  limit  as  a-*00  ,  the  thickness 

of  the  transition  layer  approaches  zero  and  the  classical 
elastic-solid-viscous'-f  luid  equations  are  recovered. 

Equations  (17)  and  (19)  can  be  re-cast  in  dimensionless 
form  and  become 

By^  ^  +  r*  R,,l  =  0  ’  (!7a) 

H(y)  =  -iR  +  (1+iR)  exp(-y/*)n  .  (19a) 


H  is  defined  to  be  equal  to  K  non-dimensionallzed  with  respect 
to  the  shear  modulus  G.  For  simplicity,  it  is  assumed  that  the 
fluid  and  the  solid  have  the  same  density;  this  is  a  good 
approximation  for  a  rubber  type  of  coating  Immersed  in  water. 

R  is  defined  to  be  the  square  of  the  ratio  of  the  speed  of 
shear  disturbances  of  the  fluid  to  that  of  the  elastic  solid.  I 
is  defined  as  the  ratio  of  nominal  coating  thickness  &  to  the 
fluid  shear  wavelength  Xf,  thus 


R  a  (cf/cs)' 
T  «  2irJl/Xf 


(20) 

(21) 
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The  fluid  shear  wave  speed  cf  and  the  elastic  solid 
transverse  speed  c8  are  given  by 


2 

cf  =  yw/p 


cs  =G/P 


The  parameters  R  and  f  can  also  be  written  as 


(22: 


(23) 


R  *  yu/G  =  (X,/Xc  )' 


r  *  a)i/c.f  =  i  (p«u/y) 


(20a) 


(21a) 


In  the  limiting  conditions  (i.e.  y«l  and  y  »l  ),  the 
unified  system  of  equations  reduce  to  the  proper  equations. 
For  y«£ ,  H  equals  to  unity  and  equation  17a  reduces  to 


+  r2Ru,  *  o  .  (2' 

For  y»i,  H  equals  -iR  and  the  differential  equation  becomes 


-J  +  ir2  u,  =  o 
,2  1 


11.4  Numerical  Solutions 


The  governing  unified  differential  equation 


&  <h  ar  >  ♦  r2R  u,  •  o 


(17a) 
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is  subject  to  the  wall  boundary  condition  and  the  imposed 
particle  velocity  at  y«h;  thus 


c 

o 

II 

o 
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(26) 

u,  (h)  =  1 

• 

(27) 

Equation  17a  subject  the  boundary  conditions  26  and  27  is 
a  split  boundary  value  problem*  The  technique  utilized  to 
solve  the  differential  equation  is  the  Runge-Kutta  numerical 
Integration  starting  at  y-0  and  going  out  to  y*h.  The  initial 
value  of  the  slope  of  ux  at  y_0  is  then  iterated  until  the 
outer  boundary  condition  ui(h)“l  is  satisfied.  The  Iteration 
technique  used  is  a  linear  interpolation  and  in  most  cases,  6-8 
Iterations  are  required  to  obtain  accuracy  of  better  than  0.1X. 
The  step  size  required  for  the  numerical  Integration  is  a 
function  of  the  parameters  R  and  T and  varies  from  y«0.01  to 
y>0.001.  In  general,  the  smaller  step  size  is  required  for 
either  thicker  coatings,  thinner  transition  layers  or  higher 
frequencies . 

It  should  be  noted  that  all  the  parameters  are  complex 

numbers,  so  that  the  iteration  on  the  initial  slope  of  u^  is 
really  a  double  iteration,  one  on  the  real  part  of  the 
derivative  and  the  other  on  the  imaginary  part. 

The  complete  copy  of  the  computer  program  used  to  solve  the 
system  of  equations  can  be  found  in  the  appendix. 


III.  DISCUSSION  AND  RESULTS 


The  results  from  the  numerical  Integration  of  the 
equations  is  presented  in  the  form  of  displacement,  shear 
stress  and  power  distribution  from  the  disturbance  location, 
through  the  transition  layer  and  in  the  elastic  coating. 

The  dimensionless  parameters  that  have  been  studied  are: 

1.  location  of  disturbance  -  h 

2*  thickness  of  transition  layer  -  n  -, 

3.  thickness  of  coating  -  &  - , 

4.  property  of  coating  -  R  -. 

The  location  of  the  disturbance  is  non-dimenslonallzed 
with  respect  to  the  nominal  coating  thickness  %  .  Thus  values 
of  h  (fig.  1)  less  than  unity  correspond  to  disturbances 
originating  within  the  coating,  while  values  of  h  greater  than 
unity  correspond  to  disturbances  within  the  fluid. 

Values  of  n  are  lnversley  related  to  the  thickness  of  the 
transition  layer.  Figure  4  shows  the  variation  of  the  real  and 
imaginary  part  of  H  (Eq.  19a)  for  two  different  values  of  n. 

A  value  of  n  equal  to  8  thus  corresponds  to  a  transition 
layer  starting  at  y/i  ■  0.75  and  ending  at  about  y/2* 1.35.  The 
transition  layer  is  thus  seen  to  be  about  equal  in  thickness  as 
the  main  coating  itself.  A  value  of  n  equal  to  32  instead 
corresponds  to  a  transition  layer  starting  at  y/1-0.95  and 
ending  ay  y /£■  1.05;  the  thickness  of  the  transition  layer  is 
thus  equal  to  10Z  of  the  compliant  coating.  The  transition 


layer  has  been  defined  as  starting  where  the  value  of  the  real 
part  of  H  is  less  than  90Z  of  its  value  within  the  compliant 
coating.  In  a  similar  manner,  the  outer  edge  of  the  transition 
layer  is  defined  to  be  the  point  where  the  value  of  the 
imaginary  part  of  H  is  within  10Z  of  the  fluid  value.  The 
point  where  the  real  part  of  H  equals  the  imaginary  part  is  the 
location  where  the  transition  layer  switches  behavior,  from 
solid-like  to  fluid- like. 

Very  rigid  coatings  with  relatively  large  shear  moduli 
will  give  rise  to  smaller  values  of  the  parameter  R  than  softer 
coatings.  The  value  of  T  ,  being  linearly  related  to  the 
dimensional  coating  thickness  "£ "  will  also  give  an  indication 
of  the  relative  magnitude  of  the  thickness  of  the  coating.  As 
it  can  be  verified  by  the  corresponding  definition  of  both  R 
and  r  >  higher  frequencies  for  the  disturbance  will  result  in  a 
simultaneous  Increase  in  both  parameters.  Table  I  gives 
typical  values  for  frequencies,-  coating  thickness  end  shear 
wave  speed  of  the  coating  for  the  rain^e  of  values  of  R  and  T 
considered  in  this  report. 
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Table  1  -  Typical  Values  of  R  and  T  for 
and  Elastic-Coating-Water  System. 


The  results  of  the  numerical  integration  of  the 
differential  equation  are  presented  in  terns  of  particle 
displacement  |u  |  ,  local  shear  stress  |t|  and  local  power 
distribution  |  P  |  for  the  different  parameters  (n,  h,  R,  T  ) 
considered.  The  power  transmitted  from  the  location  of  the 
disturbance  to  any  point  of  Interest  is  obtained  from  the 
definition 


du 


(28 


where  the  Integral  is  over  one  cycle  of  the  disturbance  and  T 
and  u  are  the  local  shear  stress  and  particle  displacement 
respectively.  Noting  that 


x  *  t (y )  exp(-lwt) 


(29 


du  *  vdt  =  -Iwudt 
T  *  1/f 


(30 

(31 


equation  28  can  be  written  as 

P  *  t (y )  u ( y )  .  (32 

Figures  5-7  show  the  particle  displacement  distribution 
for  different  location,  h,  of  the  disturbance  and  for  three 
values  of  T.  Only  the  profiles  from  the  wall  to  the  disturb¬ 
ance  is  shown.  Values  of  h  less  than  unity  correspond  to 
disturbances  within  the  solid.  The  displacement  distribution 
for  all  points  within  the  solid,  as  expected,  is  thus  seen  to 
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be  linear.  For  all  values  of  T  considered,  the  shear  fluid 
disturbance  decays  rapidly  within  the  fluid  and  then  once  it 
has  passed  the  transition  layer,  goes  to  zero  at  the  wall 
linearly.  Figure  8  shows  the  same  type  of  result  in  both  the 
linear  and  the  log  scale.  The  dotted  lines  shown  in  each  of 
the  figures  correspond  to  the  classical  solution  in  the  absence 
of  a  transition  layer.  Note  that,  while  in  the  linear  scale, 
the  difference  between  the  two  displacement  distribution  seems 
negligible,  when  the  results  are  presented  on  a  log  scale,  the 
difference  between  the  results  with  and  without  the  transition 
layer  is  more  obvious.  Figure  9  shows  the  effect  of  the 
transition  layer  thickness  on  the  displacement  distribution, 
where  it  is  seen  that  a  10Z  transition  layer  thickness  (l.e. 
n-32)  will  still  produce  a  particle  displacement  distribution 
an  order  of  magnitude  larger  than  the  classical  no  transition 
layer  solution.  Note  that  for  this  particular  case  the 
displacement  distribution  within  the  solid  is  about  10  orders 
of  magnitude  smaller  than  the  initial  value.  Figure  10  shows 
the  effect  of  the  elastic  layer  property  R  on  the  displacement 
distribution  and  is  compared  with  the  classical  solution.  As 
in  the  other  cases  the  coating  displacement  is  magnified  by  the 
transition  layer. 

Figures  11-16  show  the  shear  stress  distribution  within 
the  transition  layer  and  the  compliant  coating  for  the 
different  cases  considered.  Fig.  11  shows  the  distribution  for 


5- 


three  different  source  location.  The  distribution  for  the 
transition  layer  case  Is  equal  to  the  classical  analysis  up  to 
the  point  where  the  visco-elastic  layer  switches  from  fluidlike 
to  solid-llke;  at  this  point  the  shear  stress  suddenly 
Increases  by  an  order  of  magnitude  and  it  then  remains  constant 
to  the  wall.  The  magnitude  of  this  shift  is  dependent  upon  the 
thickness  of  the  transition  layer.  The  comparisons  between  the 
results  with  and  without  the  viscoelastic  layer  depend  upon 
whether  in  the  classical  solution  the  elastic  is  assumed  to 
exist  up  to  ymSL  or  out  to  the  location  where  the  layer  switches 
from  solid-like  to  fluid-like  (i.e.  y»1.37A  for  R-  10”5).  In 
either  case,  the  present  model  is  seen  to  cause  larger  stresses 
within  the  compliant  layer.  Figure  12  shows  the  effect  of  the 
thickness  of  the  viscoelastic  layer.  Figures  13  and  14 
illustrate  the  effect  of  different  coating  materials  for  two 
different  layer  thicknesses  and  it  is  seen  that  in  both 
instances,  the  softer  compliant  coatings  (i.e.  larger  values  of 
R)  will  result  in  greater  differences  between  the  coating-fluid 
systems  with  and  without  the  transition  layer.  Figures  IS  and 
16  show  the  effect  of  the  compliant  coating  thickness  upon  the 
shear  stress  distribution.  From  these  figures  it  is  seen  that 
the  importance  of  the  transition  layer  is  more  clear  for 
thicker  compliant  coatings  (i.e.  larger  values  of  T ) .  The 
thicker  viscoelastic  layer,  not  only  causes  the  sharp  reversal 
in  the  stress  distribution,  but  under  certain  conditions, 

causes  a  double  reversel  to  develop  (see  R-10“^*  r«  200  case). 

The  parameter  that  determines  the  magnitude  of  the 
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Interaction  between  the  fluid  and  the  coapliant  coating  is  the 
aaount  of  power  (Eq.  32)  transmitted  from  the  fluid  into  the 
coating...-  Figure  17  shows  the  effect  of  transition  layer 

4 

thickness  on  the  power  level,  and  figure  18  considers  different 
coating  materials.  In  all  cases,  the  traneltion  layer  allows 
more  power  to  be  transmitted  into  the  coating.  From  the 
results  presented  in  figure  18  it  is  seen  that  while  the 
transmitted  power  for  the  no  transition  layer  varies  linearly 
as  a  function  of  R,  the  presence  of  the  transition  layer 

reduces  the  dependence  to  about  R*^» 

Figures  19  and  20  compare  the  exponentially  varying 
transition  layer  inhomogenelty  with  a  linear  inhomogenelty  with 
the  same  end  conditions.  A  stress  driven  fluid  dlaturbance 
propagating  into  the  coating  is  shown  in  figure  21  and  22  and 
as  expected  it  is  seen  to  have  exactly  the  same  distribution  as 
the  displacement  driven  disturbance.  This  latter  result  is  as 
expected  from  the  reciprocal  theorem  of  linear  media. 

Ih.  result,  b,  P. der.ee  ,t  .l3  *r#  ln  I*"”*1  •»««■«« 
with  the  results  obtained  ln  this  report. 
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IV.  CONCLUSION 


In  this  present  analysis,  the  Interaction  between  a  fluid 
transverse  disturbance,  a  purely  elastic  compliant  coating  and 
an  Inhomogeneous  transition  layer  at  the  Interface  has  been 
studied.  The  fluid  disturbance  has  been  assumed  to  take  place 
in  the  absence  of  a  mean  flow  field,  and  either  a  transverse 
displacement  or  a  shear  stress  have  been  assumed  to  be  one 
dimensional  In  space  and  sinusoidal  In  time. 

The  results  obtained  In  this  report,  have  shown  that  a 
thin  viscoelastic  inhomogeneous  layer  can  radically  alter  the 
Interaction  between  a  fluid  shear  disturbance  and  a  compliant 
surface.  Depending  upon  the  coating  and  the  transition  layer 
thickness  as  well  as  the  frequency  of  oscillation,  the  presence 
of  the  transition  layer  lets  a  greater  amount  of  energy  Into 
the  coating. 

Although,  to  more  fully  understand  the  Interaction 
mechanism,  the  compliant  coating  was  assumed  to  be  purely 
elastic,  the  results  can  be  carried  over  to  the  case  of  a 
damped  coating.  The  transition  layer  will  allow  more  of  the 
energy  Into  the  coating  which  would  then  be  absorbed. 

The  inhomogeneity  of  the  layer  can  be  thought  as  being  due 
to  one  or  more  of  four  possibilities:  (1)  layer  with  contin¬ 
uously  varying  conditions,  (2)  a  series  of  thinner  homogeneous 
layers  each  having  different  properties,  (3)  the  coating  Is 
hydrophobic  (4)  the  coating,  being  forced  In  a  random  manner  by 


the  fluid  effectively  spreads  out  the  discontinuities  of  the 
interface  into  a  layer  with  finite  thickness. 

The  greater  percentage  of  power  transfer  from  the  fluid 
into  the  coating,  achieved  by  scans  of  a  transition  layer  aay 
stimulate  new  ideas  which  could  lead  to  more  successful 
compliant  coatings. 
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IFOUXU.C-E.l.CElSJGC  TC  451' 

*X***X*X*  ITERATION  XX  XXXXVrXXXX 


COMMENTS 


IN  CPOEP  TC  ITERATE*  ONE  MUST  HAVE  THREE  PUNS  IN 
A  GIVEN  SET.  THE  FIRST  RUN  REPRESENTS  THE  2EROTH 
ITERATION*  WHILE  THE  NEXT  TWO  ARE  PUNS  IN  WhICh 
THE  PHIAL  GUESS  FOP  DU  AT  Y  =  G  IS  VARIED  IN  SOKE 
N  ANNE c  SO  THAT  SLOPESTOP.  DERIVATIVES)*  NEECEC  IN 
CREEP  TC  CALCULATE  NEW  GUESSES*  CAN  BE  OBTAINED. 

THESE  SLOPES  AFE  OLTAINED  AS  SHOWN  IN  THE  PROGRAM  6EL0 


BETA  s' 


SMI, SMS  =  SLOPES  OBTAINED  FOR  THE  FIRST  AND 
SECOND  RUNS. 

$ N 1  * S N 2  =  SLOPES  OBTAINED  FOR  THE  FIRST  AND 
THIRD  RUNS. 

ALPHA  =  THE  VALUE  BY  WHICH  THE  INITIAL  VALUE 
OF  THE  REAL  PART  OF  DU  AT  Y  r  f  IN 
THE  FIRST  RUN  IS  TO  BE  INCREMENTED  IN 
ORDER  TO  OBTAIN  A  NEW  GUESS  FC«R  DU  AT 

Y  =  C. 

BETA  =  THE  VALUE  BY  WHICH  THE  INITIAL  VALUE 

CF  THE  IMAGINARY  PART  OF  DU  AT  Y  :  0  IN 
THE  FIRST  RUN  IS  TC  PE  INCREMENTED  IN 
CFDEF  TO  OETAIN  A  .NEW  GUESS  F0»  CU  AT 

Y  =  0. 

SLCFES  ARE  CALCULATED  ONLY  FOR  EVERY  OTHER  ITERATION. 
SLOPES  OBTAINED  FROM  THE  FIRST  THREE  RUNS  ARE  USED' 

TC  OBTAIN  THE  NEXT  TWO  ITERATIONS,  AFTER  WHICH  THE 
PROGRAM  AUTOMATICALLY  RUNS  TWO  EXTRA  RUNS  SO  THAT 
NEW  SLOPES  CAN  BE  CALCULATED.  THESE  THEN  APE  USED  TO 
OBTAIN  THE  NEXT  TWO  ITERATIONS*  AND  SC  ON.  THE  PARAMET| 
NITER  DETERMINES  THE  NUMBER  OF  ITERATIONS. 


CO  480  ►  =1*M TEF 


i 


t* 


?  1 7 
215 
215 
22  C 
221 
225 
222 
254 
221 
22k 
227 
22? 
2  25 

22  C 
5  71 
r-2  2 

223 

234 

kwC 

5  36 
2  37 
2  35 
2  79 
?*C 
5  4  1 

242 

243 

244 

245 

246 
5*7 
2*5 
2  4C 

2  c  C 
2M 
252 
?fi 

254 

s 

23  6 
2f  7 
lie 

255 
26C 
251 
262 
?f  3 

2  c  4 

265 
2*  6 
2  67 
26  £ 

265 
270 
27! 

272 

273 
57* 
273 
2  76 
2  77 
275 

27? 

250 

251 

252 
567 
264 

5*5 

266 
267 
2bk 


SLOPES  ARE  CALCULATED  FCR  EVERY  OTHER  ITERATION  AFTER  THE  2EROT 
GC.  TO !3 5 P *35 5*3* 0*355 *350 *355*350*355* 350*355) *K 
25  Cl  SF1=!UU1!2)-UU1!1))/!V1<2>-V1<1>> 

SF2  =  IUU2<2)-UU?  Cn)/(VJi2)-VlCl)J 
SNI  =  (UU1  I31-UU1  < 1> >/ (V2(3>-V?(l>>  *■ 

SN2=!UU5!3)-UU2!1) )  /  !  V2  <  3  >  -  V2  <  1  >  ) 

355  LRITE16  *  36'0  >  SMI  *SM2  *  SN1  «  SN2' 

360  FCR*'ATfl00<»-»)*/*lX*»SlQPES:_!*.46J>*J»/*10Q_C*-»)i: _ 

AA A  =  1 . 0 ”UU1 1 1 )  . 

F5P  =-UU  2(1) 

CENGM=!  SM1*SN2)-!SM2*SN1 ) 

ALPHA=C ( AAA*$N2)-!BEB*SN1) l/DENOH 
EETA=! ! SFl*E6F)-!SM2*AAA>)/CENOH 
VI C 1 )s  V 1 <1)*AL*PA 
V5(l)  =  V2(l)*»fcETA 

The  SA^r  SLOPES  ARE  US  EL  FOP  TWO  CONSECUTIVE  ITERATION. 

CC  7 C < 2  f C*270*36C«370*360*370*360f 370*3£0,37C) ,K 
37  0  N J- 3 

these  values  ape  used  fcr  the  two  extra  runs  needed  tc  calculat 

NEW  SLOPES. 

V1!5)=V1 !1 >*!2. C+PER) 

V2 (5 )=V2 ( 1) 

VI! :):V1C1) 

V2! *)=V2!1)*!1.C+PEP) 

CC  TO  350 
7  5  P  N  J  =  1 
350  CUMFY=0.0 

4  +  *  +  +  +  +  *  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  ++  *  +  +  +  + 

LC  470  U=1*NU 

IF  <hK.EC.MTER.AND.UU.E&.5)GC  TO  50C 

ALL  SYVEOLS  FFCF  THIS  PCINT  ON  HAVE  SAFE  MEANING  AS  THAT 
SFECIFTED  EARLIER. 

Ml  )=YY1 
Y! 2 )=YY2 
Y(7 > -YY3 
y<*  >=V1  !  JJ) 

Y!5)=V2(JJ) 

Z-2Z_ 

>  XXX  ITERATION  runge  -  KUTTA  xxxxx 

THIS  RUNGL-KUTTA  LOCP  IS  THE  SAME  AS  THAT  USED  PREVIOUSLY. 

CC  4  5  C  N  =  1*LIM 

IF!N.EG.l)WklTE!6*41)Y*hl*H2 

IF !N. EC. 1) WRITE! 6*42) GP*GI*RR*RI*e*Z*LlM 

Hs  H  n 

IF! A. LE.NEEG.0R.N.GT.NEND)H=H1 
CALL  DKCE!DY*Y*Z»H*W*0*5) 

CCMPY1=CPFLX (Y!2)*Y!3) ) 

UVAC1=CA5S(CGMPY1) 

CCFFY?=CFPLY(Y!4),Y!5)) 

UHAG2=C AES!C0MPY2> 

IF!P.C-T.6P.C.AN’D.Y!  1).GE.2.C)GP  TO  395 

ASSUME  THAT  AT  Y=ll  WE  ARE  OUTSIDE  THE  TRANSITION  REGION 
IF! Yll). GL. 10.0)60  TO  395 
F=! !l.-PC)*i XP!-1Y!1)**F)))*RC 
GC  Tr  357 
tot  FsP  C 

397  7AU=C0MPY2*F 

TAUP=CA6S!TAU)  *• 

kORK=TAUK*UPAGl 

IP! Y!1 ) .LE. YNUM2)GC  TC  430 

TP5  *  OR*N»£fl»il  WR  ITE<6*  300)Y(1)«Y!2)*Y(3)  *UM  AG1  * 

lTAU*TALP*WOrFfF 
CO  TC  4  f  0 

4?ClTAU*,TAU';,WCRK,f  L°’C,CR*N*EC*1>WRITE<tf  3°0>  V(  1>f  V<2>  *Y,3>  *UMAG1* 
45C  CONTINUE  * 


2fc9 

2°0 

291 

292 
29? 
294 
?c5 
296 
29  7 
256 

jce 

200 

?n 

2«2 

213 

3C* 

?rs 

3  c  6 
2C7 
3C6 


X  X  X  X  X  X  xrx~'x:  x'x’x  * '<*  X  X  X  X  X  X  X  X  X  X  X~  X  x' X 

LL'l  (J)=YC2)  *  “ .  ’  " 

LU2(d)=Y(3>. 

470  CONTINUE 

♦  *♦♦♦*.»♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

IMKK.EC.NITEF  >GO  TC  500 
460  CONTINUE  - 

XXXXXXX  END  ITERATION  X  *****  * 
GO  TO  5CC 

4ct  FCF^AT^ //tEX,  /iGNITUDES  TOC  LARGE*  ITERATE  EY  HAND.*) 
ECO  IF  (*  STf  F  .EC  .N’PUNS  >GC  TO  550 
C  ?  T  C  5 
550  £TC  f 
E»  L 


1 

2 

T 

4 

c 

6 

7 

6 

c 

10 

11 

12 

12 

1* 

15 
U 

1  7 

16 
le 

2  C 
21 
2  2 
4  - 

*4 

?5 
26 
7  7 
?fi 
pc 

3  0 
31 
22 

.’3 

le 

w 

26 

27 

26 

c6 

40 

*1 

42 

43 

L  4 

45 

46 

47 

4  t 


C  JSTUr6ANCE  IN  FLUID  * 


*  * 


■  a***- 


*  * 


SHEAR  W  A  V  l  S 
VELOCITY  -  DRIVEN 
FUNCTION  DY 

T*- 1 S  PRCGRAF  SPECIFIES  THE  DIFFERENTIAL  EQUATIONS  TO  EE 
5EE  UNI  VAC  PATH-PACN  FCP  DETAILS. 

FUNCTION  DY(Y.I) 

CC|V  YON  GR»GI.PR,RI,5 

CCNPLEX  P,U (2) ,THETA,Al*A2t A3  * D YY t X  I  X I  X  ,  G AMR  A 
DIMENSION  Y ( 5 ) 

L(1 )=C«PLX< Y(2) «Y<2)  ) 

U(2)=CNPLX(Y(4),Y(5) > 

F  =  CNFLX  CO.O.-RR ) 

XIX!X=CNPLX(0.0f 1.0) 

GAKR  A  =  CRPLX<G*?*GI) 

THE  T  A  =  <  CSGPT (XIXIX*P) ) *  G AMM  A 

ASSUME  THAT  AT  Y=10  WE  ARE  OUTSIDE  THE  TRANSITION  REGION 
IF CB.GT.f 0. C .ANC.YC1 ).GE. 2.0)60  TO  5 
IF(Y(1).GE.1G.0)G0  TO  1 
ARfc  =  Y(l  )**E 
CC  TC  2 

1  ARG=1.0F*10 

2  IF(APG.GE.1.0E+0e)GC  TO  5 

A1 *  A2  *  A  2  ARF  THE  COEFFICIENTS  OF  THE  SECOND  DIFFERENTIAL 
CYY  REPRESENTS  THE  SECOND  DIFFERENTIAL  EGUATION. 
Al=c*CR-l.)*(Y(l)**<fc-l.))*EXPC-ARG) 

A2=THET#**2 

A3=<fl.-R)*EXP<-ARG))^R 
CYYi-(Al*U(2)4>A£*U€l)»/A3 
GC  TC  6 

f  CY_Y  =  -(THETA**2  )*U<  1  )/R 
fc  GC*  TG  <10'2Cv2C«4C)«I 
10  LY’=Y<4) 

CO  TO  50 
CY=Y<5> 

CC  TC  50 
2C  CY=REAL(DYY) 
or  tc 

40  CY  =  AIYAC(CY>  ) 

50  RETURN  T  4 

t  F  D 


X  X  x| 


♦  ♦  ♦  ♦ 


X  X  X  X 


SOLVED. 
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